Empirical analysis of the worldwide maritime transportation network 
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In this paper we present an empirical study of the worldwide maritime transportation network 
(WMN) in which the nodes are ports and links are container liners connecting the ports. Using the 
different representations of network topology namely the space L and P, we study the statistical 
properties of WMN including degree distribution, degree correlations, weight distribution, strength 
distribution, average shortest path length, line length distribution and centrality measures. We 
find that WMN is a small-world network with power law behavior. Important nodes are identified 
based on different centrality measures. Through analyzing weighted cluster coefficient and weighted 
average nearest neighbors degree, we reveal the hierarchy structure and rich-club phenomenon in 
the network. 

PACS numbers: 89.75.Da,89.75.Dd,89.75.-k 



I. INTRODUCTION 

The recent few years have witnessed a great devotion 
to exploration and understanding of underlying mech- 
anism of complex systems as diverse as the Internet 
social networks [3| and biological networks [3]. As critical 
infrastructure, transportation networks are widely stud- 
ied. Examples include airline 0,ll,[i,0Ji5_ship H, bus 
[i3,[il|,[i3,[i3, subway [Ti| and railway[llllB| networks. 

Maritime transportation plays an important role in the 
world merchandize trade and economics development. 
Most of the large volume cargo between countries like 
crude oil, iron ore, grain, and lumber are carried by 
ocean vessels. According to the statistics from United 
Nations the international seaborne trade continu- 
ously increased to 7.4 billion tons in 2006 with a robust 
annual growth rate of 4.3 per cent. And over 70 per cent 
of the value of world international seaborne trade is being 
moved in containers. 

Container liners have become the primary transporta- 
tion mode in maritime transport since 1950's. Liner 
shipping means the container vessels travel along reg- 
ular routes with fixed rates according to regular sched- 
ules. At present most of the shipping companies adopt 
hub-and-spoke operating structure which consists of hub 
ports, lateral ports, main lines and branch lines, forming 
a complex container transportation network system [l^ . 

Compared with other transportation networks, the 
maritime container liner networks have some distinct fea- 
tures: (1) A great number of the routes of container liners 
are circular. Container ships call at a series of ports and 
return to the origin port without revisiting each inter- 
mediate port. It's called pendulum service in container 
transportation. While bus transport networks and rail- 
way networks are at the opposite with most of buses or 
trains running bidirectionally on routes. (2) The network 
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is directed and asymmetric due to circular routes. (3) 
Lines are divided into main fines and branch lines. Main 
lines are long haul lines which involves a set of sequential 
port calls across the oceans. Sometimes long haul lines 
call at almost 30 ports. Branch lines are short haul lines 
connecting several ports in one region to serve for main 
lines. 

We construct the worldwide maritime transportation 
network (WMN) using two different network representa- 
tions and analyze basic topological properties. Our result 
shows that the degree distribution fofiows a truncated 
power-law distribution in the space L and an exponen- 
tial decay distribution in the space P. With small aver- 
age shortest path length 2.66 and high cluster coefficient 
0.7 in the space P, we claim that WMN is a smafi world 
network. We also check the weighted network and find 
the network has hierarchy structure and "rich-club" phe- 
nomenon. Centrality measures are found to have strong 
correlations with each other. 

The rest of the paper is organized as follows: in Sec- 
tion ini we introduce the database and set up the net- 
work using two different network representations. In Sec- 
tion [Ull various topological properties are studied includ- 
ing degree distribution, degree correlations, shortest path 
length, weight distribution and strength distribution etc. 
Section Uv] discloses the hierarchy structure by studying 
the weighted and unweighted clustering and degree cor- 
relations. Centrality measures correlations and central 
nodes' geographical distribution are studied in Section 
IVl Section |Vl] gives the conclusion. 



II. CONSTRUCTION OF THE NETWORK 

We get the original data from a maritime transport 
business database named Cl-online [3| which provides 
ports and fieet statistics of 434 ship companies in the 
world. The data includes 878 sea ports and 1802 lines. 
The ports are distributed in different regions and we fist 
the number of ports in each region in Table [H 
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TABLE I: Number of sea ports by major geographic region 



Region 


No of sea ports 


Africa 


96 


Asia and Middle East 


251 


Europe 


311 


North America 


61 


Latin America 


96 


Oceania 


63 


Total 


878 




the line. The shortest path length between any two nodes 
in the space P represents the transfer time plus one from 
one node to another and thus is shorter than that in the 
space L. 

Since WMN is a directed network, we extend the con- 
cept of spaces L and P to directed networks according 
to See Fig. [TJ Line A and B are two different pen- 
dulum routes crossing at the port No. 1. (a) and (b) is 
the undirected network representation, (c) and (d) is the 
respective directed version. 

Based on the above concepts we establish the network 
under two spaces represented by asymmetrical adjacent 
matrices , and weight matrices , . The 
element of the adjacent matrix A equals to 1 if there 
is a link from node i to j or otherwise. The element 
Wij of weight matrix W is the number of container lines 
traveling from port i to port j. 

We need to define the quantities used in this weighed 
and directed network. We employ fc^(j) and k^uti"^) to 
denote in-degree, out-degree of node i in the space L, 
and k^n{i) to represent undirected degree in the space L. 
Similarly kf^{i), k^^^^{i) and k^^{i) are employed in the 
space P. Hence we have 

fc»w = E4 (1) 
fco„tW = E4 (2) 



FIG. 1: Description of the space L and the space P. (a) 
and (b) are the undirected representations in the space L 
and the space P respectively, (c) and (d) are the directed 
representations in the space L and the space P respectively. 
In the space L, a link is created between consecutive stops in 
one route. In the space P all ports that belong to the same 
route are connected. Line A (solid line) and line B (dashed 
line) are two different pendulum routes sharing one common 
node: the port No. 1. 



To construct the worldwide maritime transportation 
network, we have to introduce the concept of spaces L 
and P as presented in Fig. [TJ The idea of spaces L and 
P is first proposed in a general form in [15| and later 
widely used in the study of public bus transportation 
networks and railway networks. The space L consists of 
nodes being ports and links created between consecutive 
stops in one route. Degree k in the space L represents 
the number of directions passengers or cargoes can travel 
at a given port. The shortest path length in the space L 
is the number of stops one has to pass to travel between 
any two ports. In the space P, two arbitrary ports are 
connected if there is a container line traveling between 
both ports. Therefore degree k in the space P is the 
number of nodes which can be reached without changing 



^Lw = I](4 + 4) (3) 

which also holds for the space P. 

Strength is also divided into in-strength and out- 
strength. In the space L the in-strength of node i is 
denoted by S^^{i) and out-strength is denoted by S'^ut(*)- 
Undirected strength (total strength) is represented by 
S^^{i). They can be calculated according to the following 
equations: 

SL{^)^Y.wf. (4) 
SLt{^)=J2w^^ (5) 

which also holds for S^{i), S'^j(i), 5'^„(i) in the space 
P. 

Other quantities like clustering coefficient and average 
nearest neighbors degree also have different versions in 
directed and weighted WMN. We employ cf and cf to 
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TABLE II: Basic parameters for spaces L and P. n is the 
number of nodes and m is the number of links, (fc) is the 
average undirected degree. (C) is the average unweighted 
cluster coefficient. If) is the average shortest path length. 



Space 


n 


m 


{J^un ) 


{C) 


(0 


Space L 


878 


7955 


9.04 


0.4002 


3.60 


Space P 


878 


24967 


28.44 


0.7061 


2.66 



denote the unweighted clustering coefficient of node i in 
the space L and P respectively. Analogously k^^.i ^iid 
j are used to denote the average nearest neighbors 
degree of node i in the space L and P respectively. For 
weighted WMN we add superscript W to the above quan- 
tities and consequently they become cf" , cf^, k^^^ and 

h.WP 

III. TOPOLOGICAL PROPERTIES 
A. degree distribution and degree correlations 

First we examine the degree distributions in two 
spaces. Fig. [2] shows that in-degree, out-degree and 
undirected degree distributions in the space L all follow 
truncated power-law distributions with nearly the same 
exponents. In-degree and out-degree obey the function 
P{k) - fc-i '^ before k = 20. When fc > 20 their distri- 
bution curves bend down to the function P{k) ~ k~'^'^^. 
Unweighted degree in the space L has the same exponents 
of —1.7 and —2.95 but the critical point becomes k — 30. 
Truncated power-law degree distributions are often ob- 
served in other transportation networks like the world- 
wide air transportation network China airport net- 
work Q, U.S. airport network 0| and the Italian airport 
network It is explained in jj] that the connection cost 
prevents adding new links to large degree nodes. Analo- 
gous cost constraints also exist in the maritime transport 
network. Congestion in hub ports often makes ships wait 
outside for available berth for several days, which can 
cost ships extremely high expense. Consequently new 
links are not encouraged to connect to those busy ports. 

While in the space P three degree distributions all fol- 
low exponential distributions P{k) ^ e"''. The param- 
eters are estimated to be a = 0.0117 for in-degree, and 
a = 0.0085 for out-degree, a = 0.0086 for unweighted de- 
gree. The property that degrees obey truncated power- 
law distributions in the space L and exponential distribu- 
tions in the s pac e P is identical to public transportation 
networks [lol.llll| and railway networks Particularly 
the Indian railway network [3| has exponential degree 
distributions with the parameter 0.0085 almost the same 
with in-degree and out-degree distribution in WMN. 

Next, the relation between in-degree and out-degree 
is studied. Fig. ^c) is a plot of out-degree kout vs. 
in-degree kin- They have a positive correlation under 
two spaces. Out-degree climbs when in-degree increases. 



Evidently the in-out degree correlation is very strong. 

Finally, we want to find out the relations between de- 
grees under different spaces. Fig. [2l[d) shows undirected 
degree in the space P fc^„ vs. undirected degree in the 
space L k^^. All dots are above the diagonal, indicating 
the undirected degree in the space P is larger than in the 
space L. It is understandable because of diflerent defi- 
nition of two spaces topology. In the space P all stops 
in the same route are connected which surely increases 
degree of each node. Table [TTl lists the basic properties of 
WMN in two spaces. Average undirected degree in the 
space P is 28.44, much larger than average undirected 
degree 9.04 in the space L. 



B. Line length 

Let's denote line length, i. e. the number of stops in 
one line, as /. In Fig. [3] the probability distribution of 
line length P{1) can be approximated as a straight line in 
the semi-log picture representing an exponential decay 
distribution P{1) ~ e~"' with the parameter a = 0.13. It 
indicates there are much more short haul lines than long 
haul lines in maritime transportation. Long haul lines 
use large vessels and travel long distance from one region 
to another region while short haul lines as branch lines 
travel between several neighboring ports and provide 
cargos to main lines. For example, the line consisting 
of the following ports: Shanghai-Busan-Osaka-Nagoya- 
Tokyo-Shimizu-Los Angeles-Charleston-Norfolk-New 
York-Antwerp-Bremerhaven- Thamesport-Rotterdam- 
Le Havre-New York-Norfolk-Charleston-Colon-Los 
Angeles-Oakland- Tokyo-Osaka-Shanghai, is a typical 
long haul line connecting main ports in Asia and Europe, 
calling at ports for 24 times. 



C. Shortest path length 

The frequency distributions of shortest path lengths d 
in the spaces L and P are plotted in Fig. |4l The distri- 
bution in the space L has a wider range than in the space 
P. The average shortest path length is 3.6 in the space 
L and 2.66 in the space P (see Table ITT|) . This means 
generally in the whole world the cargo need to transfer 
for no more than 2 times to get to the destination. Com- 
pared with the network size TV — 878, the shortest path 
length is relatively small. 



D. Weight and strength distribution 

Usually traffic on the transportation network is not 
equally distributed. Some links have more traffic fiow 
than others and therefore play a more important role in 
the functioning of the whole network. Weight should be 
addressed especially in transportation networks. Here we 
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FIG. 2: (a) Degree distributions in the space L obey truncated power-law distributions with almost the same exponents. The 
turning points are at fc = 20 and fc = 30 respectively, (b) Degree distributions in the space P all follow exponential distributions, 
(c) Positive correlation between in-degree and out-degree, (d) Unweighted degree under two spaces. Degree in the space P is 
larger than in the space L. 



study four properties of weighted WMN: weight distribu- 
tion, strength distribution, in-out strength relations and 
the relations between strength and degree. The results 
are displayed in Fig. [H 

First we examine weight distributions. In Fig. [5] (a) 
two weight distribution curves are approximately straight 
declining lines before w = 40. The power-law distribu- 
tions are estimated to be P{w) ^ ty-O-S^ in the space P 
and P{w) u;-0-92 ^j^g space L. 

Next, Fig. [5] (b) shows the undirected strength distri- 
butions under two spaces both obey power-law behavior 
with the same parameter. The functions are estimated 
to be P{s) 

And we also analyze the relation between in-strength 
and out-strength in two spaces. As we can see from Fig. 
[5] (c), in-out strength relation under the space L is plot- 
ted as almost a straight line while the relation under the 
space P has a slight departure from the linear behav- 
ior at large values. But clearly they are positively 
correlated. 



Finally an important feature of weighted WMN, the re- 
lations between strength and degree, is investigated. Un- 
der two spaces the relations between undirected strength 
and undirected degree are both unlinear with the slope 
of the line approximately 1.3, which means the strength 
increase quicker than the increase of degree. This is often 
occured in the transportation networks and has its im- 
plication in the reality. It's easy for the port with many 
container lines to attract more lines to connect the port 
and thus to increase the traffic more quickly. 



IV. HIERARCHY STRUCTURE 

In this section, we explore the network structure 
of WMN through studying both the weighed and un- 
weighted versions of cluster coefficient and average near- 
est neighbors degree. Hierarchy structure and "rich club" 
phenomenon are unveiled. We conjecture this kind of 
structure is related to ship companies' optimal behavior 
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FIG. 3: Probability distribution of line length. It can be 
approximated by a straight line in semi-log plot indicating an 
exponential decay. 



to well interconnected communities (high clustering co- 
efRcient), while hubs connect many vertices that are not 
directly connected (small clustering coefficient). 

C^(k) lies above C^(fc) and the average cluster coeffi- 
cient of the network in the space P is 0.7 larger than 0.4 
in the space L. This can be explained by the fact that 
in the space P each route gives rise to a fully connected 
subgraph. With high cluster coefficient 0.7 and small 
average shortest path length 2.66 in the space P, we con- 
clude that the WMN, as expected, has the small-world 
property. 

Weighed quantities for clustering and assortativity 
measures are first proposed in i23|. Through the case 
study of WAN and SCN, [20] demonstrates that the inclu- 
sion of weight and their correlations can provide deeper 
understanding of the hierarchical organization of complex 
networks. The weighted cluster coefficient is defined as: 



1 
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FIG. 4: Frequency distributions of shortest path length under 
two spaces. The distribution in the space L has a wider range 
than in the space P. 



to minimize the transportation cost known as the hub- 
and-spoke model in transportation industry. 



which takes into account the importance of the traffic or 
interaction intensity on the local triplets. And we de- 
fine C""(fc) as the weighted cluster coefficient averaged 
over all vertices with degree k. In real weighted network 
we may have two opposite cases of the relation between 
^"{k) and C{k). If C"'(fc) > C{k) in the network, in- 
terconnected triplets are more likely formed by the edges 
with larger weights. If C""(fc) < C{k) the largest inter- 
actions or traffic is occurring on edges not belonging to 
interconnected triplets. 

In Fig. [6] we report weighted cluster coefficient under 
two spaces. Evidently the weighted cluster coefficient 
C""^(fc) and C""^(fc) are both above the corresponding 
unweig hted cluster coefficient, i. e. 0""^ {k) > C^{k), 
and C^^ik) > C^{k). This indicates some closely in- 
terconnected nodes with large degrees have the edges 
with larger weights among themselves. In other words, 
high-degree ports have the tendency to form intercon- 
nected groups with high-traffic's links, thus balancing 
the reduced clustering. This is so called "rich-club" phe- 
nomenon. 



B. Assortativity 



A. Clustering 

Cluster coefficient Ci is used to measure local cohesive- 
ness of the network in the neighborhood of the vertex. It 
indicates to what extent two individuals with a common 
friend are likely to know each other. And C{k) is de- 
fined as cluster coefficient averaged over all vertices with 
degree k. 

We plot C{k) in Fig. H Either in the space L or in 
the space P, C(fc) exhibits a highly nontrivial behavior 
with a decay curve as a function of degree A:, signaling a 
hierarchy structure in which low degrees belong generally 



There is another important quantity to probe the net- 
works' architecture: the average degree of nearest neigh- 
bors, knnik), for vertices of degree k. Average nearest 
neighbors degree of a node i is defined as: 

knn,i — ^ ^ ^ij^j (^) 
' j 

Using knn,i, one can calculate the average degree of the 
nearest neighbors of nodes with degree fc, denoted as 
knn{k). The networks are called assortative if knn{k) 
is an increasing function of k, whereas they are referred 
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FIG. 5: (a) Weight distributions in two spaces, (b) Undirected strength distributions in two spaces, (c) Relation between 
in-strength and out-strength in two spaces, (d) Correlations between strength and degree in two spaces. The slopes both equal 
to 1.3 approximately indicating nonlinear relationships in the spaces L and P. 



to as disassortative when fc„„ (fc) is a decreasing function 
oi k. As suggested in [20], weighted version of average 
degree of nearest neighbors is calculated by: 

Kn,i = -^^O^ijW^jkj (9) 

] 

From this definition we can infer that fc™„ ^ > fc„„ ^ if the 
links with the larger weights are pointing to the neighbors 
with larger degrees and ^ < knn,i in the opposite case. 

Both the weighted and unweighted average degree of 
nearest neighbors under two spaces are plotted in Fig. [7l 
The curve of fc^„(fc) lies above the curve of k^^ik). The 
k^^{k) and k^j^{k) grow with the increase of degrees at 
small degrees but decHne when degrees are large. The 
unweighted network exhibits assortative behavior in the 
small degree range but disassortative behavior in large 
range. 

When we turn to the k^^{k) and /cjf,f (fc), the weighted 
analysis provides us a different picture. We can see that 
under two spaces the weighted average degree of nearest 



neighbors exhibits a pronounced assortative behavior in 
the whole k spectrum. Since the number of WMN nodes 
is 878, this conforms with the empirical finding in 
that public transport networks are assortative when the 
number of nodes in the network N > 500 and disassor- 
tative when N < 500. 

Prom the above discussion we can see that the inclu- 
sion of weight changes the behavior of cluster coefficient 
and average degree of nearest neighbors. This property is 
identical to the worldwide airport network [23| and North 
America airport network j2l|. In both the airline trans- 
portation and maritime transportation networks, high 
traffic is associated to hubs and high-degree ports (air- 
ports) tend to form cliques with other large ports (air- 
ports). Their similar organization structure may have a 
similar underlying mechanisms. We conjecture that this 
is related to the hub-and-spoke structure which is widely 
adopted in practice by airline companies or ship com- 
panies to achieve the objective of minimizing the total 
transportation cost [H, Hi, HI • 

Fig. [8] describes a typical hub-and-spoke structure 
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FIG. 6: Cluster coefficients. C(k) under two spaces both ex- 
hibit nontrivial behavior with decay curves as functions of de- 
gree k. And weighted versions of cluster coefficient are larger 
than unweighted cluster coefficient, indicating high traffic 
edges between interconnected vertices. 




FIG. 7: Average degree of the nearest neighbors as functions 
of k. Inclusion of weight changes the behavior of knn{k): as- 
sortative behavior in the small degree range but disassortative 
behavior in large degree range, to definite assortative behavior 
of fcJJ^„ (fc) in the whole k spectrum. 



which consists of three interconnected hubs and other 
nodes allocated to a single hub. In maritime transporta- 
tion main liners travel between hubs handling large traffic 
while branch liners visit the hub's neighboring ports to 
provide cargo for the main lines. This structure allows 
the carriers to consoHdate the cargo in larger vessels to 
lower the transportation cost. This simple structure has 
the similar property of ^"(fc) > C(fc), A:™„(fc) > A:„„(fc) 
with that of WMN. And it also displays rich-club phe- 
nomenon. We think it worths investigating the relations 
between ship companies' optimal behavior and the real 
transportation network's hierarchy structure and rich- 




FIG. 8: A classical hub-and-spoke structure in maritime 
transportation. There are three central vertices having very 
strong links (high traffic) with each other and several nodes 
having weak links (low traffic) with hubs. It has the same 
property of C""{k) > C{k),k^„{k) > k„„{k) with that of 
WMN. 



V. CENTRALITY MEASURES 

In this section we analyze two centrality measures in 
social network analysis [26]: degree and betweenness. 
The most intuitive topological measure of centrality is 
given by the degree: more connected nodes are more im- 
portant. The distribution and correlations of degree has 
been discussed in Section IV. 

Betweenness centrality is defined as the proportion of 
the shortest paths between every pair of vertices that 
pass through the given vertex v towards all the shortest 
paths. It is based on the idea that a vertex is central 
if it lies between many other vertices, in the sense that 
it is traversed by many of the shortest paths connection 
couples of vertices. Hence we have 



{N - 1){N 



2) ^ 



rijk 



(10) 



where rijk is the number of shortest paths between j and 
k, and njk(i) is the number of shortest paths between j 
and k that contain node i. 

Correlations between two centrality measures are pre- 
sented in Fig. [9ljb). In both the spaces there is a 
clear tendency to a power-law relation with degree k: 
C^{k) ^ fc" with a = 1.66 in the space L and a = 1.93 
in the space P. The power-law correlations between de- 
gree and betweenness is also found in bus transportaiton 
networks [lol . [T3 | and ship transport network It's 
worth noting that this power-law relations together with 
the truncated scale-free behavior of the degree distribu- 
tion implies that betweenness distribution should follow 
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FIG. 9: (a) Probability distribution of betweenness obeys 
truncated power-law distribution. (b) Betweenness is a 
straight line in log-log picture indicating a power-law rela- 
tion with degree. The exponent in the space L is estimated 
to be 1.66 and that in the space P is estimated to be 1.93. 



a truncated power law. This behavior is clearly identi- 
fied in Fig. [DJa). We find the betweenness centrality has 
two-regime power-law behavior P(C^) ^ C~". For the 
two spaces, exponents are almost the same: a = 0.14 at 
small degree regime and a = 1.0 at large degree regime. 

The power-law relations between degree and between- 
ness suggest that they are consistent with each other. It 
is proved in the comparison of each port's degree and 
betweenness. The 25 most connected ports are hsted in 
Table. IIIII Singapore is the most busy ports in the world 
with the largest degree and betweenness. Antwerp and 
Bushan are the second and third either in degree or in 
betweenness measures. Only 5 ports in these ports are 
not listed in the 25 most central ports in betweenness 
measure. WMN is not like the case of the worldwide air- 
line network Q which has anomalous centrality due to 
its multicommunity structure. The difference may due 
to the fact that there are less geographical and polit- 
ical constraints in maritime transportation than in air 



TABLE III: The 25 most connected ports in the worldwide 
maritime transportation network. * These ports are not 
among the 25 most central ports. 



Rank Ports 


Degree 


Betweenness 


Region 


1 


Singpore 


120 


124110.1258 


Asia 


2 


Antwerp 


102 


113368.6161 


Europe 


3 


Bushan 


92 


69094.7490 


Asia 


4 


Rotterdam 


87 


78097.8754 


Europe 


5 


Port Klang 


83 


62111.6226 


Asia 


6 


Hongkong 


78 


46072.9799 


Asia 


7 


Shanghai 


75 


37748.4316 


Asia 


8 


Hamburg 


60 


40362.4625 


Europe 


9 


Valencia 


60 


27346.0956 


Europe 


10 


Le Havre 


58 


48231.1636 


Europe 


11 


Gioia Tauro 


55 


20148.4667 


Europe 


12 


Yokohama 


54 


19716.7287 


Asia 


13 


Kaohsiung 


52 


15363.5781 


Asia 


14 


Port Said* 


52 


16121.8476 


Africa 


15 


Bremerhaven 


50 


24998.9088 


Europe 


16 


Colombo* 


48 


15128.362 


Asia 


17 


Tanjung Pelepas 


4fi 






18 


Jeddah* 


45 


8570.7724 


Middle East 


19 


Jebel Ali 


44 


20137.5255 


Middle East 


20 


Ningbo* 


44 


12315.9515 


Asia 


21 


Algeciras 


43 


18701.2619 


Europe 


22 


Barcelona* 


43 


15256.3903 


Europe 


23 


Kobe 


43 


8869.7313 


Asia 


24 


New York 


42 


25935.6692 


North America 


25 


Kingston 


41 


23281.3223 


Latin America 



transportation. Ships can travel longer distance than 
airplanes. And airports are usually classified into inter- 
national and domestic airports and international airlines 
are hmited to connect international airports instead of 
domestic airports. So there are distinct geographically 
constrained communities in WAN. In the maritime trans- 
portation there are no such constraints. Sea ports basi- 
cally are all international ports with the possibility to 
connect to any other sea ports in the world. 

In Fig. [To] we plot the 25 most connected ports on 
the world map. They show unbalanced geographical dis- 
tribution mainly located in Asia and Europe, includ- 
ing 13 ports in Asia and Middle East, 1 in Africa, 9 
in Europe, 1 in North America and 1 in Latin Amer- 
ica. Particularly they are located along the east-west 
lines. Lines in maritime transportation are usuahy di- 
vided into east-west lines, north-south lines and south- 
south lines [l§ . [27| . The fact that 25 most connected 
ports in the world are in east- west trade routes represents 
rapid growth and large trade volume in Europe- America, 
Asia- America and Asia-Europe trade [iTj . 
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VI. CONCLUSION 

In this paper we have presented an empirical study of 
the worldwide maritime transportation network (WMN) 
under different representations of network topology. We 
study the statistical properties of WMN and find that 
WMN is a small world network with power law be- 
havior. There are strong correlations in degree-degree, 
strength-degree and betweenness-degree relations. Cen- 
tral nodes are identified based on different centrality mea- 
sures. Based on the analysis of weighted cluster coeffi- 
cient and weighted average nearest neighbors degree, we 
find that WMN has the same hierarchy structure and 



"rich-club" phenomenon with WAN. We conjecture that 
this structure is related to optimal behavior both exist- 
ing in air transportation and maritime transportation. 
So our future research direction is the evolution model- 
ing of WMN using optimal behavior to reproduce real 
properties in WMN. 
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